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Using tlie non-linear Ginzburg-Landau (GL) tlieory, we obtain the possible vortex configurations 
in superconducting thin films containing a square lattice of antidots. The equilibrium structural 
phase diagram is constructed which gives the different ground-state vortex configurations as function 
of the size and periodicity of the antidots for a given effective GL parameter k* . Giant-vortex states, 
combination of giant- and multi- vortex states, as well as symmetry imposed vortex-antivortex states 
are found to be the ground state for particular geometrical parameters of the sample. The antidot 
occupation number Uo is calculated as a function of related parameters and comparison with existing 
expressions for the saturation number Us and with experimental results is given. For a small radius of 
antidots a triangular vortex lattice is obtained, where some of the vortices are pinned by the antidots 
and some of them are located between them. Transition between the square pinned and triangular 
vortex lattices is given for different values of the applied field. The enhanced critical current at 
integer and rational matching fields is found, where the level of enhancement at given magnetic 
field directly depends on the vortex-occupation number of the antidots. For certain parameters of 
the antidot lattice and/or temperature the critical current is found to be larger for higher magnetic 
fields. Superconducting/normal H ~ T phase boundary exhibits different regimes as antidots are 
made larger, and we transit from a plain superconducting film to a thin-wire superconducting 
network. Presented results are in good agreement with available experiments and suggest possible 
new experiments. 
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I. INTRODUCTION 

Superconducting (SC) samples with periodic arrays of 
pinning sites have received much attention over the last 
decade. It is now well established that these artificial 
pinning centers (i) hold great potential for enhancing the 
critical parameters of the sample and (ii) give rise to 
different kinds of vortex behavior that is not observed in 
the presence of random pinning. In this respect, arrays of 
microholes (antidots)i'^'^''^'tiSiSiL2i2ii2iiiiiS and submicron 
magnetic dotS j-^^d'^d^ have been studied, as their pres- 
ence in the SC film strongly modifies the vortex structure 
compared to the one in non-patterned films. ^^•-'^^ 

Direct imaging experiments,^ magnetization 
and transport measurements and theoreti- 

cal simulationsi2ii2iS2i2iiS^ of vortex structures in 
samples with periodic pinning centers have shown 
that the vortices form highly ordered configurations 
at integer ff„ = n^o/S and at some fractional 
Hp^q — ^^q/S (n,p,q being integers) matching fields, 

where $o = hc/2e = 2.07 • lO~'^Gcm^ is the flux 
quantum, and S is the area of the primitive cell of the 
artificial lattice. This remarkable variety of stabilized 
vortex lattices may even be broadened by multiple 
possible degeneracies. These commensurability effects 
between the pinning array and the vortex lattice are 
responsible for an enhanced pinning and consequently 
increased critical currents. Very recently Karapetrov 
et alm^ investigated vortex configurations in a single 
crystal superconducting heterostructures with an array 



of submicron normal metal islands by scanning tunneling 
microscopy. They observed the coexistence of strongly 
interacting multiquanta vortex lattice with interstitial 
Abrikosov vortices. Different vortex phase transitions 
are given, which occur when the number of magnetic 
fiux quanta in the sample changes. 

Motivated by those experimental studies on perforated 
superconductors, significant efforts have been made on 
the theoretical side as well. For example, extensive 
molecular dynamics simulationsiSii^iSSiSiiS^ in the Lon- 
don limit have been performed in an attempt to calcu- 
late the vortex structure and their dynamics in a peri- 
odic pinning potential. Although the general behavior of 
vortex lattices was accurately described, made approx- 
imations are valid only in certain range of parameters. 
Namely, in the London approach, vortices are consid- 
ered as classical point-particles (with different models for 
their interaction) and the influence of the antidots is in- 
troduced through model hole potential, which in princi- 
ple can never be generalized. Recently, Nordborg and 
Vinokur^^ discussed in the detail interaction of vortices 
with an arbitrarily large cavity, but still within the Lon- 
don theory. This study was actually an extension of 
the work of Mkrtchyan and Shmidt^S^ who crudely esti- 
mated the maximum possible number of vortices trapped 
by a single insulating inclusion with radius i? as = 
R/2^{T), where £,{T) is the temperature dependent co- 
herent length. For regular arrays of pinning centers the 
saturation number becomes Ug ~ (i?/^(0))^ due to the 
vortex- vortex interactions The antidot- vortex interac- 
tion and the following maximal occupation number of 
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each antidot appear to be crucial for many phenomena. 
For example, experiments on thin films with a lattice of 
holes showed a "localization transition"i^: all vortices 
drop inside the holes when the coherence length becomes 
larger than the interhole spacing. In Rcf.^ it was shown 
that the antidot size realizing the optimum pinning is ac- 
tually field-dependent. The effective vortex-pinning po- 
tential and saturation number of the pinning sites for 
different temperature and applied dc fields were recently 
investigated experimentally by means of ac-susceptibility 
measurements, for superconducting films with an array of 
antidots^'* and for the case of not fully perforated holes 
{i.e. blind holes). ^'^ 

Most of the experiments on perforated superconduct- 
ing films are carried out in the effective type-II limit 
(k* = A^/d^ ^ I/a/S, d being the thickness of the super- 
conducting film and A the magnetic penetration depth). 
In this regime, the vortices act like charged point par- 
ticles and their interaction with periodic pinning poten- 
tial can be described using molecular dynamic simula- 
|-jQj^g^^8ji9j20j2i^^ However, the overlap of vortex cores 
(with size ^) , and the exact shape of the inter- vortex 
interaction (depending on the superconducting material 
properties reflected through k), may significantly modify 
the vortex structures and consequently the critical cur- 
rent when this criteria is no longer satisfied. 

Besides, the vortex-pinning and the critical current en- 
hancement, higher critical field (Hcs) near an open cir- 
cular hole in a thin film (the so-called "surface supercon- 
ductivity") has been predicted theoreticalljiS^ and con- 
firmed experimentally.^^ Cusps in the H — T boundary 
were observed, which occurs when the number of vortices 
which nucleate inside the hole increases by one, similarly 
to the known Little-Parks effect. The ratio between the 
critical fields in perforated samples was estimated in lim- 
iting cases: Hc3/Hc2 — 1 when R ^ (or R << ^) and 
Ha,/ Hc2 — 1.695 when R-^ oo. 

In this work superconducting films with square ar- 
rays of antidots are treated within the phenomenolog- 
ical Ginzburg-Landau theory. This approach considers 
vortices as extended objects and no approximations have 
to be made on e.g. the vortex-vortex interaction and/or 
the vortex-antidot interaction. In Sec. II, the details of 
our numerical formalism are given. Sec. Ill deals with 
vortex lattices in perforated films in homogeneous mag- 
netic field, with emphasis on the number of pinned and 
interstitial vortices as function of the antidot-size and in- 
terhole distance. In case of weak pinning potentials, i.e. 
small size anti-dots, we discuss the triangle to square vor- 
tex lattice transition in Sec. IV. In Sec. V, we address 
the behavior of critical current in the sample as function 
of the applied field, for different geometrical parameters, 
and temperature. The dependence of the critical field on 
temperature, and different regimes in the H — T phase 
diagram are discussed in Sec. VI for different antidot- 
size. All presented findings are then summarized in Sec. 
VII. 



II. THEORETICAL FORMALISM 

In this work, we consider a thin superconducting film 
(of thickness d <C ^, A) with a square array of holes (ra- 
dius R, period W) immersed in an insulating media in 
the presence of a perpendicular uniform applied field H 
(see Fig. QJ. To describe the superconducting state of 
the sample we solve the coupled nonlinear GL equations, 
which are written in dimensionless units in the following 
form22i^: 

(~i\' - ly ^ = ^ {I - \^\^) , (1) 



--K*Ai'= — («'*V* - W**) - l^-l^ A (2) 

We measure the distance in units of the coherence length 
^, the vector potential A in c?i/2e^, the magnetic field H 
in Hc2 = ch/2e£^^ = k\^Hc, and the order parameter 
in \J —aj fi with a, /3 being the GL coefficients. 

The magnitude of the applied magnetic field H — 
n$o I S va determined by the number n of fiux quanta 
$0 — hc/2e = 2.07- 10~^Gcm^ piercing through the rect- 
angular simulation area S = Wx x Wy = N^NyW"^, with 
Wx(y) = Nx{y)W and and Ny are integers. At the su- 
perconductor/insulator interface we impose the bound- 
ary condition corresponding to zero normal component 
of the superconducting current. The periodic bound- 
ary conditions for A and ^' (simulating the periodicity 
of both superconducting film and antidot lattice) have 
the forin2^ 

A{p + h)=A{p)+^^l^{p), (3) 



*(/9 + hi) = * • exp{2^^ir]^{p)/<^o), (4) 

where bi (i = x,y) are lattice vectors, and rji is the 
gauge potential. These boundary conditions imply that 
A, ^' are invariant under lattice translations combined 
with specific gauge transformations ry^.y. Other quanti- 
ties, such as the magnetic field, the current or the or- 
der parameter density are periodic. We use the Landau 
gauge ^0 = HxCy for the external vector potential and 
r]x = HWxU + Cx, Vy = Cy, with C^, Cy being con- 
stants.^^ Without antidot lattice, when the film is invari- 
ant under infinitely small translations, the free energy 
does not depend on Cx, Cy. The vortex lattice is only 
shifted relative to the simulation region when Cx and Cy 
are varied. This is not the case for an antidot lattice, 
when the change of these parameters leads to a displace- 
ment of the vortex lattice relative to the holes, leading 
to a variation of the free energy. In general, one has to 
minimize the free energy with respect to Cx, Cy. It can 
be shown that such a minimization gives a zero current 
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FIG. 1; Schematic view of the sample: a superconducting 
film (thickness d) with a regular array (period W) of circular 
antidots (radius R). 



when averaged over the cell area. We find that for a su- 
percell having one hole, the optimal values are given by 



C. 



0, ±7r and that for the supercell with 2 holes the 
choice Cx,y = provides the minimum free energy. 

We solved the system of Eqs. H1I2|I self-consistently 
using the link variable approach^ in a finite-difference 
representation of the order parameter and the vector 
potential using a uniform cartesian space grid {x,y). 
The first GL equation is solved with a Gauss-Seidel it- 
eration procedure. '^'^ The vector potential is then ob- 
tained with the fast Fourier transform technique. The 
temperature is indirectly included in the calculation 
through the temperature dependence of the coherence 
length ^{T) ^ ^{0)/ ~ T/Tao\ and penetration depth 
A(r) = A(0)/^|1 -T/Tcol, where T^o is the critical tem- 
perature at zero magnetic field. 
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III. VORTEX LATTICES - INFLUENCE OF 
GEOMETRICAL PARAMETERS 

We first consider a supercell containing four holes 
(^Wx = Wy = 2W) of radius R, with lattice period W 
(see Fig. ^ . Although our approach is valid for any inte- 
ger number of fiux-quanta piercing through the simula- 
tion region, we will restrict ourselves here to the so-called 
(integer and fractional) matching vorticities. 

Figure 12 shows contour plots of j^'p of the vortex lat- 
tice in case of antidots of radius R — ^ and interhole 
distance W — 8.0^ for different matching fields The 
film thickness is d = 0.1^ and the effective GL parameter 
K* — }? /dS^ = 10. At the first matching field all vortices 
are trapped in the antidots [Fig. [2Ja)]. Because of their 
small radius, each antidot is able to pin only one vortex, 
and additional vortices localize at interstitial sites when 
H > Hi [see Fig. EJb) for H = H3/2]. At the second 
matching field [Fig. ^c)] vortices occupy all interstitial 
sites, forming again a square lattice. For H = [Fig. 
|2Id)] vortices form an ordered lattice with an additional 
vortex at every other interstice, added to the H2 case. 
Note that the size of the vortex cores at neighboring in- 
terstitial sites differs: two vortices at the same interstice 
strongly interact and effectively bound each others core- 
areas, while the neighboring single interstitial vortex does 
not suffer from any lack of space resulting in its larger 



FIG. 2: Contour plot of the Cooper-pair density (white/dark 
color - low/high density) in the superconducting film with 
antidots of radius R — ^ and interhole distance W = 8.0^ for 
the matching fields: //„ = 1 (a), 3/2 (b), 2 (c), 5/2 (d), 3 (e), 
4 (f), 9/2 (g), and 5 (h). The film thickness is d = 0.1^ and 
the effective GL parameter k* = 10. Dashed lines indicate 
the antidot lattice. 

core. At the third matching field (three vortices per unit 
cell) two interstitial vortices in adjacent cells alternate in 
position [see Fig. [2Ie)], preserving the two- fold symme- 
try, but the vortex unit cell is 4 times the antidot lattice 
unit cell. At H — H4 [Fig. Hff)] we observe the first 
evidence of the competition between the Abrikosov vor- 
tex lattice (characteristic for thin film superconductors) 
and the symmetry of the pinning lattice, as vortices (in- 
cluding pinned ones) form a slightly deformed hexagonal 
lattice. Notice that for H = Hq^2 [Fig-I21g)] the number 
of vortices per antidot unit cell is n = 9/2 and the vortex 
lattice unit cell is twice the antidot unit cell. At the fifth 
matching field [Fig. Efh)] the dense packing of vortices 
and their consequent strong interaction with the anti- 
dot lattice result in the restoration of the square lattice 
symmetry but the vortex lattice is tilted over 35° with 
respect to the antidot lattice. Our results are in excel- 
lent agreement with the experiment of Refi and previous 
molecular dynamics simulations,— in certain parameter 
range. Namely, the vortex configurations are mainly de- 



FIG. 3: Contour plot of the Cooper-pair density (left column) 
and phase of the order parameter (right column, limited to a 
single antidot lattice unit) for the sixth matching field H = 
and for different values of the hole radius: R = 0.8^ (a), 1^ 
(b), 1.6C (c), 1.7C (d) and 2.1^ (e). The lattice period is 
W = and the GL parameter k* = 10. 



termined by the pinning force of each antidot, and the 
vortex-vortex interaction. The latter is very dependent 
on the density of vortex packing, as the known expres- 
sions for the vortex-vortex interaction do not take into 
account possible overlap of vortex cores. On the other 
hand, the antidot-vortex pinning potential is determined 
by the antidot-size and the period of the antidot-lattice. 
Our approach takes all these aspects into account and 
their influence will be discussed in the remainder of this 
paper. 

As shown in Fig. |3 for small hole radii, each antidot 
pins only one vortex (the hole occupation number rio, the 
number of vortices sitting in the holes, equals 1), and the 
remaining vortices reside between the holes. One expects 
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FIG. 4: Contour plot of the Cooper-pair density (left column) 
and phase of the order parameter (right column) for the fourth 
matching field H — H4 and for the different values of the 
lattice period: W = 4.2^ (a), 4.4^ (b), 5^ (c), 5.2^ (d) and 5.8^ 
(e). The radius of the holes is R = 0.8^ and GL parameter 
k' = 10. 



that, for larger hole radius R, vortex configurations with 
multi-quanta vortices in each hole can become energet- 
ically preferable Figure 13 shows the contour plots 
of the Cooper-pair density at the sixth matching field 
for different antidot radii. The antidot lattice period is 
W = 6^, the thickness is d = O.lf and n* = 10. The 
number of vortices captured by each hole changes from 
one for R = 0.8^ [Fig. E^a)] to five for R = 2.1^ [Fig. 
|3Je)]. The vortex arrangement outside the holes is de- 
termined not only by their mutual interaction, but also 
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FIG. 5: Contour plots of the Cooper-pair density for k* = 10 
(a-c) and k* = 0.1 (d-f). The lattice period is W = 6(, and 
the radius is i? = 1^. Figures (g-i) show the phase of the 
order parameter of the states shown in (d-f). The first row is 
for H = H2, the second for H — H3 and the bottom row for 
H = H4. 



by the attraction with the antidots and the repulsion by 
their pinned vortices. For small radius R a multivortex 
structure is found at the interstitial sites, as apparent 
from Fig.|3{a,b) (see also Fig.|21). By further increasing 
R some of the vortices enter the holes and the remain- 
ing vortices arc strongly caged between the antidots, re- 
sulting in the formation of giant vortices [Fig. El^d)] and 
symmetry imposed vortex-antivortex pairs^^ [Fig. EIc)]. 
This is apparent from the contour plot of the phase of 
the order parameter 4* (right column of Fig. |3J . 

Similar behavior can be achieved if the hole-size is kept 
the same, but the period of the hole lattice is decreased, 
as illustrated in Fig. 0] For small period W the intersti- 
tial vortices form a giant vortices [Fig. ^a)] because of 
the strong interaction with the pinned vortices in the 
antidots. With increasing W one extra vortex is de- 
pinned and, due to the symmetry of the sample, vortex- 
antivortex pair with four vortices and one antivortex is 
formed in each interstitial site [Fig. Efb)]. Further in- 
crease of W leads to a triangular vortex structure at the 
interstitial sides with chosen orientation that minimizes 
the energy between neighboring cells [Figs. ^c-e)]. 

It is well known that the vortex-vortex interaction 
changes sing at the point k = 1/ \/2. For k > 1/ -\/2, vor- 
tices repel each other while for k < 1/V2 they attract. To 
see how this attractive interaction modifies the different 
vortex lattice configurations we consider a sample with 
small K. Fig. [S] shows the contour plots of the Cooper- 
pair density for k* = 10 (type-II regime) and k* = 0.1 
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FIG. 6: Contour plots of the Cooper-pair density (left col- 
umn) and the phase of the order parameter (right column) of 
the sample in Fig. |^ for different metastable vortex states at 
H = Hi and k* = 0.1. 



(type-I regime) for the second, third and fourth match- 
ing fields. For the given parameters of the sample and 
for K* = 10 each hole pins one vortex and the remaining 
vortices sit at interstitial sites [Figs. |H a-c)]. The occupa- 
tion number of each hole is increased to two in the type-I 
sample [see Figs. EJd-i)]) due to the enhanced expulsion 
of the magnetic field by the superconductor. Moreover, 
because of the attractive interaction between vortices, 
giant vortices become energetically more favorable con- 
trary to the case for n* = 10. Due to the instabilities of 
vortex states, which is common for type-I superconduc- 
tors, variety of metastable vortex structures can be found 
(see also Refi^). As an example we show in Fig. differ- 
ent metastable vortex states of the sample in Fig. for 
K* — 0.1. The free energies of those states are: F/Fq —- 
0.3268 (a), -0.2823 (b), -0.2787 (c) and -0.2755 (d). The 
ground state free energy [Fig. [SIf,i)] is F/Fq = —0.3759. 
Notice that because of the attractive interaction a giant 
vortex state is always favored. 
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FIG. 7: The vortex phase diagram: the dependence of the 
hole occupation number Uo on the radius of the holes R and 
the distance between them W, for k* = 10 (solid curves) and 
K* =0.1 (dashed curves) at the matching fields H = H2 (a), 
H = ff4 (b) and H = //g (c). The gray region shows the 
normal state region for k* = 10 (tC =0.1 when delimited by 
the full curve). 
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FIG. 8: (a) The antidot occupation number as a function of 
the antidot radius from the London theory— Us = R/2^{T) 
(solid line) and from the GL theory for the period W = 5^ 
(dashed curve), W = 10^ (dotted curve), W = 15^ (dash- 
dotted curve) and W — 30^ (gray curve) for the fourth match- 
ing field and k* = 10. (b) The critical hole radius correspond- 
ing to Uo = 1 ^ no = 2 transition as a function of the period 
W . Solid curve is obtained from the London theory^- and 
dashed (dotted) curve is the result from the GL theory for 
H ^H2 {H = H4) and k* = 10. 



To summarise the above findings, we constructed the 
equilibrium vortex phase diagram, which shows the de- 
pendence of antidot-occupation number Uo on R and W 
for two values of k* , at the second, fourth and sixth 
matching fields (Fig. EI). The ground- and metastable 
states are determined in our calculation by comparing 
the energy of all stable vortex states found when start- 
ing from different randomly generated initial conditions. 
The procedure of finding the free energy of the differ- 
ent metastable states was similar to that used for the 
case of mesoscopic superconducting disksi^Si^ It should 
be noted, that an energetically unfavorable state remain 
stable in the wide range of variation of R and W. There- 
fore, the transitions between the vortex states with differ- 
ent occupation numbers are of first-order. It is seen from 
this figure that Uo increases as the applied field is in- 



creased, which is in agreement with experimental results 
and theoretical predictionsi2iiSiS4S& With decreasing k* 
the threshold hole-radius for capturing another vortex de- 
creases due to the smaller pinned vortex-interstitial vor- 
tex repulsion. 

Let us compare our numerical results for the hole- 
occupation number rio with existing theoretical predic- 
tions. The saturation number Hs {us = Uq for larger 
fields) is usually estimated as Ug = R/2(,{T)m^ Fig.lSfa) 
shows the hole occupation number Hq obtained from this 
expression and the one from our GL calculation for dif- 
ferent period of the antidots. It is seen from this figure 
that this expression underestimates no for small period 
W (dashed and dotted curves). This is due to the fact 
that the last expression does not account for the interac- 
tion between vortices sitting at different holes. For larger 
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period W the occupation number is smaller in our calcu- 
lation for a given radius of the holes (gray curve). A more 
accurate analysis was presented by BuzdinS£ for bulk su- 
perconductors within the London approach. However, 
his estimation of the critical hole radius R/^ « (W/^)'^^'^ 



(for < A) and R/^ ; 



,2/3 



(for A <C W) corresponding 



to the transition from single flux-quantum to two flux- 
quanta captured by the hole, differ from our numerically 
exact results, i.e., the magnitude of the critical hole ra- 
dius is largely overestimated in RefiS£ for both small and 
large period W [see Fig. |HIb)]. 

The maximum number of flux quanta that can be 
trapped in a pinning center in a thin superconducting 
film was recently studied experimentally using scanning 
Hall probe microscopy^i and ac susceptibility measure- 
ments!^ In the latter case the saturation number was ob- 
tained from the transition to different dynamic regimes, 
as the interstitial vortices have higher mobility than those 
pinned by the antidots. They studied thin Pb films con- 
taining a square antidot array of period d = 1.5/im. The 
antidots had circular (square) shape with radius R — 
330nm (size a — 0.8^m), the film thickness was d — 80nm 
{d — lOOnm) and the coherence length at zero tempera- 
ture was estimated ^(0) = 30nm (^(0) = 33nm) in Ref>ii 
(RefiS). Let us first discuss the results for the sample 
of Ref;^, where the experimentally obtained saturation 
number was 71^=2 at T/Tcq — 0.77. Fig.l^a) shows the 
antidot occupation number function of temper- 

ature for different applied matching fields. At small ap- 
plied fields {H < H3) the occupation number is equal 
to two, which is in agreement with the experimentally 
obtained Ug. With increasing applied field H > one 
more vortex is trapped by the holes, i.e. Uo = 3, which is 
now larger than the experimental value. At higher tem- 
peratures T > 0.89Tco, rio again becomes equal to two. In 
this case one would estimate the saturation number from 
n, w i?/2C(T)25 to heus^l for 0.967Tco < T < T^o and 
Us = 2 for O.SeSTco <T < 0.967Tco. We found the occu- 
pation number equal to Uq — I only for the second match- 
ing field at the temperature range 0.985Tco < T < Tc. 
The estimation of Buzdin'^^ for the critical hole radius 
R^ < ^(r)A(T)^, where the transition from tIq = 1 to 
Uo — 2 occurs, gives the temperature range T < 0.985Tco. 
We found this transition at this temperature only for the 
second matching field. For larger fields the occupation 
number is always larger than unity. The giant vortex 
state is found only at H = H4 for T > 0.984rcO and 
the vortex-antivortex state is formed at H = H5 for the 
temperatures T > 0.986Tco. 

Up to now we use the temperature dependence for the 
coherence length and penetration depth as shown at the 
end of Sec. II which is obtained from the BCS theory^ 
and is valid near T^. In this case the GL parameter k 
is temperature independent. In recent experiments on 
Pb arrays of nanowires arrays Stenuit et alm^ found that 
the following temperature dependence of the coherence 
length C(r) = ^{0)V\'^-t^\/{l - t^) and penetration 
depth A(T) = / , which leads to a tempera- 
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FIG. 9: The antidot occupation number Uo as a function of 
temperature T /Tco for different matching fields for circular 
antidots {R = 0.33/xm) (a) and square (a = 0.8/im) (b) holes. 
The lattice period for both samples \s W = 1.5//m and film 
thickness is ci = 80nm (a) and d = lOOnm (b). The inset in (a) 
shows the hole occupation number Uo of the same sample as 
a function of temperat ure for a different temperature depen- 
dence of C(r) = ^(0)v'|l -t*|/(l-t^) and k = K(0)/(l + t^), 
where t = T/Tco and k(0) = A(0)/C(0). 



ture dependence of the GL parameter k = k(0)/(1 -I- 1^), 
agrees better with experiment. Here t = T/Tco and 
k(0) = A(0)/^(0). These expressions are obtained from 
the two-fluid model. We calculated the hole occupation 
number Uo using the above temperature dependencies, 
which is shown in the inset of Fig. Ela). It is seen from 
this figure that the transition from Uo ~ n to Uo = n-\ 
occurs now at higher temperatures, but the results are 
qualitatively similar with the earlier results. 

Fig. Elb) shows the occupation number func- 
tion of temperature for different matching fields for the 
sample of RefiS. The experimentally obtained saturation 
number was ris = 3 for temperatures T/Tcq > 0.97. Our 
calculations give the same occupation number Hq — 3 for 
this range of temperatures but only for fields H < H4. 
The expression for the saturation number Hs ~ a/4^(T) 
gives in this case Ug — I. At larger fields the occupation 
number increases, but still there will be interstitial vor- 
tices in the sample»22i These interstitial vortices lead to 
a larger dissipation in the sample which was used as the 
criterium for the determination of the saturation num- 
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FIG. 10: The ground-state vortex lattice at _ff = ^^9/2 as 
a function of the antidot-radius R and their periodicity W, 
for ft* = 10 (a) and k* = 0.1 (b). The solid lines denote 
the first order transitions between the states with different 
antidot-occupation number Uo, and dashed ones depict sec- 
ond order configurational transitions. The insets show the 
Cooper-pair density plots of the corresponding states indi- 
cated by the numbers in the phase diagram. 



ber Us- But our calculations show that the appearance 
of interstitial vortices does not indicate the saturation of 
trapped vortices. 

We have shown in our recent paper-"* that a rich va- 
riety of ordered vortex structures: a combination of gi- 
ant vortices with mltivortices and vortex-antivortex pairs 
are found in perforated superconducting samples for frac- 
tional matching fields. Here we consider the dependence 
of these vortex states on the effective GL parameter k* . 
As an example we constructed the equilibrium vortex 
phase diagram for H = Hg/2 rational matching field (i.e. 
4.5 flux quanta per antidot) as a function of R and W, 



and for k* = 10 [Fig. HHIa)] and k* = 0.1 [Fig. UUJb)]. In 
type-II regime [Fig. llOf a)]. for larger period W the vor- 
tex configuration always consists of individual vortices, 
except for Uo flux lines pinned by each hole. One ex- 
tra vortex per two holes is shared between the adjacent 
cells (insets 2, 5) or situated in every other cell (inset 
7). With decreasing W, the interstitial vortices become 
strongly caged between the neighboring antidots, result- 
ing in the distortion of the individual-vortex lattice. In 
this case, the n-th matching field becomes larger than the 
second critical field for W/^ < \/2nn. But the supercon- 
ducting state in perforated films still survives due to en- 
hanced superconductivity in close proximity around the 
holes (due to surface superconductivity) . If the radius of 
the holes is then increased, individual vortices captured 
at interstitial sites can merge for H > Hc2 to a giant 
vortex [(region 1 in Fig. Iior a)]. This transition does not 
show any hysteretic behavior and is, therefore, of second 
order (similarly to the case for mesoscopic disks^S). The 
creation of these giant vortex states is favored because of 
the repulsion of the vortices by the supercurrents around 
the holes (and consequent compression of vortices in the 
central part of the interstitial regions) . 

The influence of caging depends on the number of 
confined vortices; the combinations of giant- and multi- 
vortices may be formed in the interstitial sites (insets 1, 
4 and 6). For interstitial vorticity 3, a vortex-antivortex 
pair may nucleate, so that local vortex structure con- 
forms with the square symmetry of the pinning lattice 
[see inset (3) in Fig.HHIa) and Ref.^-*]. 

Fig. ITUT b) shows the ground-state phase diagram 
found for k* =0.1. Compared to Fig. llOf a). the thresh- 
old antidot-radius for capturing another vortex decreases 
due to the enhanced screening of the applied field. Due 
to the attractive vortex- vortex interaction in type-I sam- 
ples, giant-vortex states become energetically favorable 
at the interstitial sites and spread over the majority of 
the W — R phase diagram (light gray areas) . For a dense 
antidot lattice, giant vortices with different vorticity are 
found in adjacent cells [L = 3 and L — 4 (inset 1), and 
L = 2 and L — 3 (inset 3)]. Contrary to the type-II case, 
these giant vortices can split to smaller giant vortices for 
larger spacing of antidots. They exhibit single-vortex be- 
havior, forming the lattice of 2-quanta and single-quanta 
vortices (insets 2 and 4). Such new quasi- Abrikosov 
lattices of giant-vortices result from the competition of 
vortex-vortex attraction and imposed square symmetry 
of pinning. At the same time, these competing inter- 
actions cause the complete disappearance of the vortex- 
antivortex structures as found in type-II samples. 



IV. TRIANGULAR TO SQUARE VORTEX 
LATTICE TRANSITION IN THE PRESENCE OF 
A SQUARE ANTIDOT LATTICE 

It is well known that the regular triangular vortex lat- 
tice has the lowest energy in superconductors with no 
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FIG. 11: The phase diagram: square (white region) and de- 
formed triangular (light gray region) vortex lattice as a func- 
tion of the radius R and period W of antidots at H — ffi/g(a) 
and H = i?3/ie(b) for k* = oo (solid curves) and k* = 0.9 
(dashed curves). The insets show the Cooper-pair density 
plots of the corresponding states. 

pinningiii As we have showrn above the square lattice 
of pinning sites impose its own symmetry on the vor- 
tex structure. If the vortex-pinning strength in a peri- 
odic square array is reduced, the vortex-vortex repulsion 
starts to dominate over the pinning force and the tri- 
angular lattice is recovered. Transition between these 
phases was recently studied in Refi^ as a function of the 
amplitude of the vortex-pinning site interaction and the 
characteristic length scale of this interaction within the 
London theory (i.e. k* = oo). They showed that the 
transition between triangular and square vortex lattice 
occurs for increasing strength of the pinning potential in 
the case of small values of pinning potential length scale. 
In Ref.'^^ a model periodic pinning potential was intro- 
duced, the parameters of which are difficult to relate to 
any growth parameters of the sample. 

To circumvent the latter problem we studied the case 
of weak pinning potential by introducing small antidots. 
Thus in order to decrease the pinning force in our calcu- 
lations we just reduced the radius of the antidots R for a 
given period of antidot lattice W. Calculations are done 
for a 8 X 8 unit cell {Wx = Wy = 8W) with grid points 




x/'q x/i; xJt 



FIG. 12: Contour plot of the Cooper-pair density for different 
vortex states at, H — Hi. The radius of the antidots is R = 
0.14^, the period is W = 5^ and effective GL parameter is 

K* = oo. 



256 X 256. Fig. ^2 shows the phase diagram: the tran- 
sition between the pinned (white region) and triangular 
(light gray region) vortex lattice as a function of the ra- 
dius R and period W of antidots for the applied fields 
H = Hi/g{a) and H = _ff3/ig(b). Let us first discuss 
the results for H = Hi/^. When the pinning strength 
is small (small i?), it is energetically favorable to form a 
triangular lattice, where vortices are located between the 
holes [left inset in Fig. Illf a)]. For larger radius of the 
holes a square vortex lattice becomes the ground state 
[right inset in Fig. Illf a)] . The critical radius of the holes 
R to pin the vortices decreases with increasing period, 
contrary to the one corresponding to the case of two-flux 
quanta captured by the holes [see Fig. |HIb)] . If we de- 
crease the GL parameter k the transition between pinned 
and triangular vortex lattices decreases (dashed curves in 
Fig. 111!) due to the short range interaction between the 
vortices. 

In Ref,^^ the phase diagram for the transition between 
the triangular and square vortex lattice was found to be 
the same for all submatching fields not exceeding H — 
Hi, except H — Hi and H — i?i/2. Contrary to this 
results our calculations give a different phase diagram for 
different fractional matching fields. As an example, we 
show in Fig. Illf b) the transition lines between triangular 
and square pinned vortex configurations for H = H^/iq. 
The triangular vortex lattice is formed where some of the 
vortices are pinned and some of them are located between 
the antidots [see the inset of Fig. lTTT b)]. For this value of 
the field the triangular vortex lattice is found for larger 
values of the period W and radius R of the antidots. 

For the fields H = Hi and H — H1/2 the pinned vortex 
lattice has a square symmetry. An intermediate vortex 
configuration for these fields was obtained in Refi^ where 
vortices in odd rows of pinning centers are depinned and 
are located between the pinning sites forming a kind of 
triangular lattice. Our calculations show that such vor- 
tex configurations can be found only as a metastable state 
and for relatively larger values of the period {W > 4^). 
To make this more clear we plot in Fig. E| the ground 
state (a) and metastable (b,c) vortex configurations at 
the first matching field. In addition to the triangular 
vortex state given in Ref.^^ [see Fig. I12f b)] another tri- 
angular vortex state is found, where all the vortices are 
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located in the interstitial region [Fig. I12r c')]. 



V. CRITICAL CURRENT OF PATTERNED SC 
FILMS 

In the previous sections we showed that vortex config- 
urations that are commensurate with the periodic arrays 
of antidots exhibits well-defined matching phenomena, 
which leads to pronounced peaks in the critical current 
(see for example Ref4i). However, the stability of these 
vortex states strongly depend on the parameters of the 
sample. For example, a multi-quanta vortex state be- 
come energetically favorable for large radius of the holes, 
while small holes can capture only a single vortex. The 
additional vortices located in interstitial sites reduces the 
critical current considerably. Therefore, we first investi- 
gate the critical current of our sample as a function of 
the relevant antidot parameters. 

The first step to calculate the critical current is to ac- 
curately determine the vortex ground state for given ap- 
plied magnetic field, in a manner described in previous 
section. Then the applied current in the x direction is 
simulated by adding a constant Acx to the existing vector 
potential of the applied external field. With increasing 
Acx we find a critical value of A^x such that a stationary 
solution to Eqs. (|1I2II cannot be found since a number 
of vortices is driven in motion by the Lorentz force. The 
current jx in the sample corresponding to the given value 
of Acx is obtained after integration of the x-component 
of the induced supercurrents in the y-cross-section. The 
maximal achievable value of jx denotes the critical cur- 
rent jc- 

Fig. 1131 shows the critical current density jc (in units 
of jo = ci/c2^/47rA^) as a function of applied magnetic 
field (normalized to the first matching field Hi) for dif- 
ferent values of the antidot radius R for fixed value of 
the antidot lattice period W . For small radius (solid cir- 
cles), where only one vortex can be pinned by the hole, 
the peaks at the matching fields decrease with increasing 
applied field. The opposite behavior is found when there 
is a caging effect?^ i.e. jc{Hn) < jc{Hn+i), which e.g. is 
found for radius R = 1.9^ (open circles) for Ug — 1 and 
Uo — 2. This effect occurs when there are the same num- 
ber of interstitial vortices but the number of pinned ones 
are different at the different matching fields. In this case 
the interstitial vortices feel a stronger repulsive interac- 
tion when there are a larger number of pinned vortices. 
As is shown in Fig. E| (open circles), a higher critical cur- 
rent is found for the third matching field, when a double 
vortex occupies each hole and a single one is located at 
the interstitial, than for the second matching, with one 
vortex in each hole and a single interstitial vortex (see 
the insets of Fig. ll3|l . This effect disappears with further 
increasing the radius R due to the different occupation 
number Uo, i.e. no interstitial vortices at H ~ H2. 

In order to show the range of radius R and period 
W of antidots, where this caging effect is active, we 
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FIG. 13: Critical current density (in units of jo = 
C-ffc2C/47rA^) as a function of the applied magnetic field (in 
units of the first matching field Hi) for three values of the 
antidot radius: R — (sohd circles), R = 1.9^ (open cir- 
cles) and R — 2.5^ (squares). The antidot lattice period is 
W = 10^, the film thickness is d = 0.1^ and the effective 
GL parameter is k* — 10. The insets show the contour plots 
of the Cooper-pair density at the second (1) and third (2) 
matching fields for R = 1.9^. 



constructed a R — W phase diagram for H — H2 and 
H = H3, shown in Fig. I14r a'). The shadowed area indi- 
cates the vortex state with a single interstitial vortex for 
both H = H2 (solid line) and H = H-^ (dashed curves). 
Fig-EJb) shows the ratio jc{Hz) / jc{H2) as a function of 
period W . The critical radius R is taken from the middle 
of the region (dotted curve). It is seen from this figure 
that, although we have the same vortex structure for all 
values of the period 4^ < < lOf , the enhancement 
of jc is found only for W ^ 6.6f . For small period the 
pinned vortices at H = H3 suppresses superconductivity 
around the holes and interstitial vortices are easily set 
into motion, reducing the critical current. 

Fig. 1 151 shows the critical current density as a function 
of the field for two values of the period: 1^ = 4^ (open 
circles) and W = 8^ (solid circles) at i? = 1.25^. As we 
showed above, the jc{H) curve shows pronounced max- 
ima at integer fields Hi, H2 and H^ and at some of the 
fractional matching fields. However, while the qualitative 
behavior oi jc{H) in Fig. E|is as expected, its quanti- 
tative behavior reveals a counterintuitive phenomenon. 
Namely, one expects higher critical current in the sample 
with larger interhole distance, simply due to the pres- 
ence of more superconducting material. Indeed, that is 
the case for H < Hi, where the superconductor is able to 
compress all flux lines in the holes. However, for higher 
magnetic fields, the critical current drops sharply imme- 
diately after the first matching field Hi, which is related 
to the appearance of interstitial vortices. On the other 
hand, the smaller interhole distance affect the hole occu- 
pation number (see Fig. [7|l, and the additional vortices 
after H = Hi are still captured by the holes (as illus- 
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FIG. 14: (a) The dependence of the antidot occupation num- 
ber rio as a function of antidot radius R and lattice period 
W (see also Fig. |7|l. The solid (dashed) line indicates the 
transition between the states with different Uo at H = H2 
{H = Hi). The shadowed area indicates the vortex state with 
a single interstitial vortex for both applied fields (the occupa- 
tion number in this region is rio = 1 for H = H2 and rio = 2 
for H = H3). (b) The ratio of jc{Hz) / jc{H2) as a function of 
the period W , where the radius of the antidots correspond to 
the dotted line in the middle of the dashed area. 



trated by Cooper-pair density plots in the inset of Fig. 
I15|l . Consequently, the critical current in this case is 
larger for smaller periodicity. Note that even for smaller 
periodicity a sharp drop in jc is observed for H > Hi, 
as every additional vortex disturbs the stability of the 
vortex lattice. Even at H = H2, although all vortices are 
captured by the holes, the critical current is lower, due to 
a stronger suppression of the order parameter around the 
holes compared to the H ~ Hi case. The height of the 
matching peaks is decreasing with further increasing field 
(due to the presence of interstitial vortices), which agrees 
with experiment (see Ref.^), and these peaks strongly di- 
minish for higher fields as the vortex-flow overwhelms the 
pinning potential. 

When we apply a dc current into the superconductor 
the vortex lattice is distorted before the vortices start 
moving. To illustrate this phenomenon, we plot in Fig. 1161 
the Cooper pair density of the superconducting film at 
the applied currents (in y direction) j = (first column), 
j = 0.5jc (second column) j — 0.95jc (third column) for 
different matching fields. At the first [Fig. llGf a)] and 
second [Fig. [TBT b')] matching fields all the vortices are 
displaced over the same distance, conserving the square 
symmetry in the lattice of vortices. At larger fields, when 
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FIG. 15: Critical current density of the superconducting film 
as a function of the applied magnetic field for two values of 
the antidot lattice period: W = 8^ (solid circles) and W = 
(open circles). The insets show contour plots of the Cooper- 
pair density at the second (1,3) and third (2,4) matching fields 
for W = AS, (1,2) and 1^ = 8^ (3,4). The radius of the holes 
is i? = 1.25^, the film thickness is d = 0.1^ and the effective 
GL parameter is k* = 10. 



there is a large number of interstitial vortices [Fig. I16f c- 
e)], the vortex configuration is changed by the current 
and some of the vortices are jammed at the interstitial 
sites. If we initially have giant vortices [Fig. Ilfcir fl] they 
can be split into multivortices with increasing j. Our 
calculations also show that there is no transition from 
the multivortex state to the giant vortex state when we 
increase the applied current, and the occupation number 
of the antidots Uq is found to be independent of j. 

Another interesting feature following from the dis- 
placement is found for fractional matching fields. For 
example. Fig. [2td) shows alternating two- vortex - single 
vortex structure at H = Hc^ /2 1 where applying small cur- 
rent in y-direction can shift the excess-vortex from one 
interstitial site to another. Note that resulting state has 
identical configuration and energy as the previous one. 
In order to estimate the energy barrier between these 
two vortex states we performed calculations for a super- 
conducting film of thickness d — 13nm with an array 
of antidots with period W = 1/im, radius R = 0.13/im, 
at temperature T = 0.9Tc. We take ^(0) = 40nm and 
A(0) = 80nm, which are typical values for Pb thin films. 
We found an energy barrier of Ai^ = 6.2meV, which 
is significantly higher than the thermal activation en- 
ergy at this temperature {kT = 0.56meV), but still low 
enough for successful switching by a relatively weak cur- 
rent. Moreover, when an ac current is applied to the 
sample, the vortex can shift back and forth between the 
adjacent cells, resulting in resonant dissipation. 

So far, we presented results at a fixed temperature. In 
what follows, we include temperature in our numerical 
analysis through the temperature dependence of the co- 
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FIG. 16: Contour plot of the Cooper-pair density for H — 
i/i (a), H = H2 (h), H = Hs (c,g), H = H4 (d,f), and 
H = H5 (e), and for the applied current j = (first column), 
j — 0.5 jc (second column) and j = 0.95jc (third column). 
GL parameter is k* = 10, the period of the antidots lattice is 
W = 6^ and the radius of the antidots is J? = 0.8^ (e), R = 1^ 
(a-d), R = 1.3C (f) and = 2^ (g). 



herence length ^. We now consider the superconduct- 
ing film with thickness d = 20nm, interhole distance 
W = 1/xm, and antidot radius R = 0.2^m. We choose 
the coherence length ^(0) = 40nm and the penetration 
depth A(0) = 42nm, which are typical values for Pb films. 
Fig.llTlshows the calculated critical current of the sample 
as a function of the applied field normalized to the first 
matching field at temperatures T/Tco — 0.86 0.98. As 
expected, decreasing the temperature leads to a larger 
critical current for all values of the applied field. The 
relative height of the peak at zero field with respect to 
one at the first matching field increases with increasing 
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FIG. 17: Critical current density of the perforated super- 
conducting film as a function of the applied magnetic field 
(in units of the first matching field H\) at temperatures 
T/Tco = 0.86 ^ 0.98. The lattice period is W = l/xm, the 
antidot radius is _R = 0.2/im, and film thickness is d = 20nm. 



temperature [see Fig. ITM a)]). At higher temperatures, 
i.e. for ^(T) > i?, a certain suppression of the order 
parameter is present around the antidots [right inset in 
Fig. I18r a)] as the core of pinned vortices overlaps with 
the interstitial regions. Consequently, the suppressed or- 
der parameter leads to a smaller jc- The caging effect is 
found for temperatures T < 0.93Tco [see Fig. I18r b)] and 
it disappears with temperature when approaching T^q, 
since the vortices entirely cover the interstitial regions 
and effectively destroy superconductivity. 

This effect that the critical current is larger for larger 
fields was recently observed experimentallyi^ The consid- 
ered sample was a Pb film of thickness 50nm, with square 
antidots of size a — 0.5/im and period W = 1.5/im. The 
coherence length and the penetration depth at zero tem- 
perature were estimated to be ^0 = 40nm and Aq = 
80nm. Although plotted for other purposes, Fig. 6(b) 
in Refj-S demonstrates a clear overshoot of the critical 
current aX H = H3 with respect to the one at H = H2, 
at the temperature 0.974Tco. Fig. 1191 shows the compar- 
ison of the calculated critical current density (dots) with 
experiment (solid line). Our jc{H) curve shows the same 
qualitative behavior as the experimental one, though a 
quantitative agreement is lacking for the experimentally 
estimated values of ^0 and Aq. Better correspondence was 
achieved for smaller values of indicating somewhat 
"dirty" sample in the experiment. No further attempts 
were made to improve the quantitative agreement with 
experiment because of the different determination of jc 
in the experiment and in our theory. In our calculations 
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FIG. 18: The mtiooiMHo) /MHi) (a) iind j,{H3)/jc{H2) 
(b) as a function of temperature for the sample with parame- 
ters given in Fig. 1171 The insets show the Cooper-pair density 
plots at temperatures indicated by the solid circles in the main 
figure (a). 



we use a dynamical criterium, i.e. we assume normal 
state as soon as vortices are set in motion, whereas in 
transport measurements a certain value of the threshold 
voltage was used to determine the critical current and 
the surface barrier at the edges is important. Therefore, 
our result should be considered as a lower limit to the 
experimental critical current. The qualitative behavior 
of jc at the matching fields should not be influenced by 
these facts. 



VI. SUPERCONDUCTING/NORMAL Tc{H) 
PHASE BOUNDARY 

The presence of antidot lattice in a superconducting 
film not only enhances the vortex-pinning, which was dis- 
cussed in the previous section, but also affects substan- 
tially the nucleation of superconductivity. Due to the su- 
perconducting/vacuum interface at the antidots, surface 
superconductivity will be important around each antidot, 
at fields above the bulk critical field Hc2 (T) . This makes 
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FIG. 19: Numerically obtained jc{H) characteristics of the 
superconducting film with antidot arrays (parameters given 
in the figure) for the coherence length at zero temperature 
^0 = 40nm (solid circles) and = 25nm (open circles). The 
solid curve denotes experimental data [taken from Ref^]. 



it possible to enhance the critical field in patterned super- 
conducting films above Hc2{T) and even beyond the third 
critical field H^^iiT). The ratio HcsiT) / Hc2{T) tends to 
the value 1.69, the enhancement factor for a semi-infinite 
slab.^^ However, for a dense antidot lattice a much larger 
enhancement can be achieved. Namely, if the antidots are 
sufficiently closely spaced, almost the entire sample may 
become superconducting at high fields through surface 
superconductivity. 

The critical field of superconducting Pb films with a 
square array of antidots was investigated in Refi^ by the 
magneto resistance measurements. The experimentally 
obtained H — T phase boundary shows a cusp-like behav- 
ior with cusps at integer and some fractional matching 
fields. The amplitude of the cusps depend on the resistive 
criterion: the cusps become sharper and their amplitude 
increases with decreasing this criterion. 

We investigated numerically the H—T phase boundary 
for a superconducting film of thickness d = 50/xm in the 
presence of a regular array of square antidots with lattice 
period W = 1.5^m. We take the coherence length at zero 
temperature as ^(0) = 40nm and penetration depth as 
A(0) = 42nm. Fig. EHl shows the calculated Tc{H) phase 
diagram for different sizes of the antidots: a = 0.25/im 
(solid curve), a — 0.5/im (dashed curve), a — 0.75/im 
(dotted curve), and a — 1.0/im (dash-dotted curve). For 
comparative reasons, we plotted also the phase boundary 
for a plain film (thin solid curve) with the same coher- 
ence length ^(0) = 40nm, obtained from the well-known 
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FIG. 20: H — T phase boundary for the superconducting film 
with an antidot array. The film thickness is d = 50nm, the 
period of the antidot lattice is = 1.5/xm, the antidot size 
is varied as a = 0.25/im (solid curve), a = 0.5/im (dashed 
curve), a = Q.75/im (dotted curve), and a = I.O/im (dash- 
dotted curve). Thin solid curve denotes the upper critical field 
[Hc2) of the plain superconducting film [Eq. J^], thin dashed 
curve gives the third critical field = 1.69Hc2 for a plain 
superconductor-vacuum boundary, and thin dotted curve is 
the critical field of a superconducting strip with thickness lo — 
0.5^m [Eq. ©]. 

expression for the upper critical field 

It can be easily seen that the antidot lattice has a pro- 
found influence on the critical magnetic field, as com- 
pared to a reference non-patterned film. The critical 
temperature is enhanced at every field, and vice versa, re- 
gardless of the size of the antidots. Note also that match- 
ing features are present in Tc{H) at integer matching 
fields. For small radius of the antidots matching peaks at 
higher integer matching fields H > H2 are weakly pro- 
nounced, due to the small hole-saturation number (see 
Sec. III). We did not observe clear evidence of fractional 
matching features. 

For small radius of the antidots the sample basically 
acts as a non-patterned film for temperatures close to Tco 
and the dependence of the critical temperature on the 
applied field is almost linear. For larger sizes of the anti- 
dots (e.g. dotted and dash-dotted curves in Fig. l^ . the 
critical field becomes substantially higher than the third 
critical field of a semi- infinite slab (thin dashed curve), 
and the peaks at matching fields are more pronounced. 
In addition, Tc{H) exhibits a parabolic background as 
for a thin slab in a perpendicular field, as well as to a 



thin film in a parallel field, which can be described in the 
London limit by^ 

(jj 

where lo stands for the width of the superconducting 
strip. 

VII. CONCLUSIONS 

We have studied the vortex structure of a thin super- 
conducting film with a regular array of antidots, which 
shows a rich variety of ordered vortex lattice configura- 
tions for different matching and fractional matching fields 
Hn- For small radius of the holes, the vortex configura- 
tions with one vortex captured in each hole and the others 
located in the interstitial sites are realized, where inter- 
stitial vortices form regular patterns, either as multi- or 
giant vortices, or combination of giant- and multi- vortex 
states. For particular geometrical parameters of the sam- 
ple and the applied field, a symmetry imposed vortex- 
antivortex configuration is found. Depending on the ratio 
between the hole radius R and the interhole distance W, 
multi-quanta vortices may be forced into the antidots, 
in spite of their low saturation number at smaller mag- 
netic fields. To illustrate the transition between possible 
multi-quanta states in the holes we showed a diagram 
of the occupation number rio as a function of the radius 
of the holes and interhole distance for different values of 
the effective GL parameter. Uo increases with decreasing 
K* due to the enhanced expulsion of the magnetic field 
from the superconductor and giant vortices become en- 
ergetically favorable because of the attractive interaction 
between the vortices. 

When the pinning force of the antidots is small, i.e. 
small radius of the antidots, the triangular vortex lat- 
tice becomes energetically favorable. Depending on the 
applied field all the vortices can be located between the 
antidots, or some of them are pinned by the antidots and 
some of them are located between the pinning centers. 
We calculated the phase diagram which shows the tran- 
sition between the triangular and pinned square vortex 
lattices for two values of the applied field and GL pa- 
rameter K. We found that the results from the simple 
London theory for the phase diagram are different from 
our GL results for different applied fields. Moreover, we 
could not find triangular vortex structures for the fields 
H = Hi and H = H1/2 as a ground state. 

The critical current jc of the sample shows well defined 
peaks at different matching Hn and fractional matching 
fields, indicating that vortices are strongly pinned by an- 
tidots. However, the level of jc enhancement at particular 
magnetic field strongly depends on the antidot occupa- 
tion number 7io. For certain parameters of the sample, 
the critical current becomes larger at higher matching 
fields, contrary to conventional behavior. 
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We also studied the Tc{H) phase boundary of regularly 
perforated superconducting film. When an antidot array 
is present the critical temperature Tc{H) is enhanced 
compared to a non-patterned film and distinct cusps 
in the phase boundary are found for different matching 
fields, which is in agreement with the experimental4i 
This behavior is in contrast to the Little-Parks^ like 
structures found in finite size superconductors 4^ The in- 
crease of the antidot size for given lattice period leads 
to the change of the Tc{H) background from linear to 



parabolic behavior except for T near Tcq. 



ACKNOWLEDGEMENTS 



This work was supported by the Flemish Science Foun- 
dation (FWO-Vl), the Belgian Science Pohcy (lUAP), 
the JSPS-ESF NES program, and the ESF-AQDJJ net- 
work. 



* Electronic address: golibjon.berdiyorov(S)ua.ac.be( 
t Electronic address: francois.peeters@ua.ac.be 

^ K. Harada, O. Kamimura, H. Kasai, T. Matsuda, A. Tono- 
mura, and V.V Moshchalkov, Science 274, 1167 (1996). 

^ M. Baert, V.V. Metlushko, R. Jonckheere, V.V. 
Moshchalkov, and Y. Bruynseraede, Phys. Rev. Lett. 74, 
3269 (1995). 

^ V.V. Moshchalkov, M. Baert, V.V. Metlushko, E. Rosseel, 
M.J. Van Bael, K. Temst, R. Jonckheere, and Y. Bruynser- 
aede, Phys. Rev. B 54, 7385 (1996). 

* V.V. Moshchalkov, M. Baert, V.V. Metlushko, E. Rosseel, 
M.J. Van Bael, K. Temst, Y. Bruynseraede, and R. Jonck- 
heere, Phys. Rev. B 57, 3615 (1998). 

5 V. Metlushko, U. Welp, G.W. Crabtree, R. Osgood, S.D. 
Bader, L.E. DeLong, Zhao Zhang, S.R.J. Brueck, B. Ilic, 
K. Chung, and P. J. Hesketh, Phys. Rev. B 60, R12 585 
(1999). 

® L. Van Look, B.Y. Zhu, R. Jonckheere, B.R. Zhao, Z.X. 
Zhao, and V.V. Moshchalkov, Phys. Rev. B 66, 214511 
(2002). 

A.V. Silhanek, S. Raedts, M. Lange, and V.V. 

Moshchalkov, Phys. Rev. B 67, 064502 (2003). 
® A.V. Silhanek, S. Raedts, M.J. Van Bael, and V.V. 

Moshchalkov, Phys. Rev. B 70, 054515 (2004). 
^ A.V. Silhanek, L. Van Look, R. Jonckheere, B.Y. Zhu, S. 

Raedts, and V.V. Moshchalkov, Phys. Rev. B 72, 014507 

(2005). 

1° S. Raedts, A.V. Silhanek, M.J. Van Bael, and V.V. 

Moshchalkov, Phys. Rev. B 70, 024509 (2004). 
" A.N. Grigorenko, G.D. Howells, S.J. Bending, J. Bekaert, 

M,J, Van Bael, L. Van Look, V.V. Moshchalkov, Y. 

Bruynseraede, G. Borghs, LI. Kaya, R.A. Stradling. Phys. 

Rev. B 63, 052504 (2001). 

M.M. Doria and G.F. Zebende, Phys. Rev. B 66, 064519 
(2002). 

" J.I. Martin, M. Velez, J. Nogues, and Ivan K. SchuUer, 
Phys. Rev. Lett. 79, 1929 (1997). 

J.I. Martin, M. Velez, A. Hoffmann, I.K. Schuller, and J.L. 
Vicent, Phys. Rev. Lett. 83, 1022 (1999). 
M.V. Milosevic and P.M. Peeters, Phys. Rev. Lett. 93, 
267006 (2004). 

A.A. Abrikosov, Zh. Eksp. Teor. Fiz. 32, 1442 (1957). 
W.H. Kleiner, L.M. Roth, and S.H. Antler, Phys. Rev. A 
133, 1226 (1964). 

C. Reichhardt, J. Groth, C.J. Olson, Stuart B. Field, and 
Franco Nori, Phys. Rev. B 54, 16108 (1997). 
C. Reichhardt, C.J. Olson, and Franco Nori, Phys. Rev. B 
57, 7937 (1998). 



C. Reichhardt, G.T. Zimanyi, and N. Gronbech- Jensen, 
Phys. Rev. B 64, 014501 (2001). 

C. Reichhardt and N. Gronbech- Jensen, Phys. Rev. B 63, 
054510 (2001). 

B.Y. Zhu, F. Marchesoni, V.V. Moshchalkov, and F. Nori, 
Phys. Rev. B 68, 014514 (2003). 

G. Karapetrov, J. Fedor, M. lavarone, D. Rosenmann, and 
W.K. Kwok, Phys. Rev. Lett. 95, 167002 (2005). 

H. Nordborg and V.M. Vinokur, Phys. Rev. B 62, 12408 
(2000). 

25 G.S. Mkrtchyan and V.V. Shmidt, Sov. Phys. JETP 34, 
195 (1972). 

'^^ M. Doria, S.C. de Andrade, and E. Sardella, Physica C 

341-348, 1199 (2000). 
2^ A. Bezryadin and B. Pannetier, Physica Scripta. 66, 225 

(1996). 

2* A.I. Buzdin, Phys. Rev. B 47, 11416 (1993). 
2® A. Bezryadin and B. Pannetier, J. Low Temp. Phys. 98, 
251 (1995). 

^° V.A. Schweigert, P.M. Peeters, and P.S. Deo, Phys. Rev. 
Lett. 81, 2783 (1998). 

V.A. Schweigert and P.M. Peeters, Phys. Rev. B. 57, 13817 
(1998). 

^2 M.M. Doria, J.E. Gubernatis, and D. Rainer, Phys. Rev. 
B 39, 9573 (1989). 

R. Kato, Y. Enomoto, and S. Maekawa, Phys. Rev. B 47, 
8016 (1993). 

G.R. Berdiyorov, M.V. Milosevic, and P.M. Peeters, Phys. 
Rev. Lett. 96, 207001 (2006). 
^5 G.R. Berdiyorov, M.V. Milosevic, and P.M. Peeters, Phys- 
ica C 437-438, 25 (2006). 

M. Tinkham, Introduction to Superconductivity (New 
York: McGraw-Hill) (1996). 

G. Stenuit, S. Michotte, J. Govaerts, and L. Piraux, Su- 
percond. Sci. Technol. 18, 174 (2005). 
P.G. de Gennes, Superconductivity of Metals and Alloys 
(Benjamin, New York, 1966). 

W.V. Pogosov, A.L. Rakhmanov, and V.V. Moshchalkov, 

Phys. Rev. B 67, 014532 (2003). 
""^ G.R. Berdiyorov, M.V. Milosevic, and P.M. Peeters, Eu- 

rophysics Letters 74(3), 493 (2006). 
"1 E. Rosseel, T. Puig, M. Baert, M.J. Van Bael, V.V. 

Moshchalkov, and Y. Bruynseraede, Physica C 282-287, 

1567 (1997). 

W.A. Little and R. Parks, Phys. Rev. A 113, 97 (1964). 
G.R. Berdiyorov, B.J. Baelus, M.V. Milosevic, and P.M. 
Peeters, Phys. Rev. B 68, 174521 (2003). 



